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DETERMINING PROPERTIES OF THE ANTENNAE SYSTEM 
- MERGING ABILITY FOR RESTRICTED N-BODY 

> , 

^ . Hanns P. Petsch^ and Christian Theis^ 



Abstract. Motivated by the closest major merger, the Antennae Galax- 
ies (NGC 4038/4039), we want to improve our genetic algorithm based 
modeling code MiNGA fTheis [l999p . The aim is to reveal the major 
interaction and galaxy parameters, e.g. orbital information and halo 
I ■ properties of such an equal mass merger system. Together with the 

O ' sophisticated search strategy of MiNGA, one needs fast and reliable 

models in order to investigate the high dimensional parameter space 
^ ' of this problem. Therefore we use a restricted JV-body code which is 

based on the approach by Toomre & Toomre (|1972|) . however with some 
refinements like consistent orbits of extended dark matter halos. Re- 
cently also dynamical friction was included to this code (Petsch [20071). 
^ \ While a good description for dynamical friction was found for mass ra- 

. tios up to q = 1/3 (Petsch & Theis [2008)l . major merger systems were 

only imperfectly remodeled. Here we show recent improvements for a 
major merger system by including mass-loss and using NFW halos. 
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1 Introduction 



I> . The Antennae System (NGC 4038/4039) is the closest observable major merger. 
vTji \ It serves as a prototype for interacting galaxies and tidally induced activity. Cur- 
. rent high precision observations over a large range of wavelengths give insight to 
' the dynamical and morphological properties. If one wants to perform state-of-the- 
art high resolution self-consistent models, a reliable set of orbital and structural 
parameters is required. Due to the high dimensional parameter space of this prob- 
lem, one needs fast and reliable models as well as a sophisticated strategy to pick 
out the best one. Therefore we use the program MiNGA (Theis IT999p , where an 
improved restricted JV-body code is coupled to a genetic algorithm. It was shown 
in Theis & Kohle (|200ip . that it is possible to reproduce the large scale features of 
an observed system with this method. A similar approach was successfully used 
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by Wahde & Donncr ((200T|) and Ruzicka et a;. (|2007|) . 



2 Method 



The introduction of dynamical friction in our restricted JV-body code allowed us 
to model close interactions and late stages of merging. Up to a mass ratio of 
q = 1/3, the orbital decay and the basic tidal features of a satellite merging with 
an isothermal halo can be reliably reproduced (Petsch & Theis I2008p . While 
a mass and distance dependent description for the Coulomb logarithm together 
with an orientation correction of the frictional force was applicable for these minor 
mergers, it did not perfectly work for a major merging system. 

In this previous work we have focused on the determination and parameter- 
isation of the Coulomb logarithm In A. We found the best description of the 
dynamical friction using a modified Chandrasekhar approach (Eq. l2.ip . 

= -FivM, a)Cf^ {VM cos if3) + sin (/?)) In A (2.1) 



dt vjj 



In A = In 



^ _l_ Mhalo(?'Af) 



M 



(2.2) 



The acceleration dvM/dt of a massive particle M ("satellite") depends on the 
background mass-density p, the mass of the perturber M and its velocity vm ■ For 
more details on the function F{vm, c), refer to Binney & Tremaine (1987) Eq. (7- 
18). Cf is a scaling factor of order unity. As galaxies have density gradients, the 
force might point not exactly opposite to the velocity. Therefore we introduced an 
orthogonal component which is adjustable via an angle (3. The Coulomb logarithm 
In A is the relation between the maximum impact parameter femax and the impact 
parameter Bq that leads to a 90° degree deflection. We found its best parameteri- 
sation to be mass and distance dependent - cf. Eq. (|2.2p . where Mhaio(?'M) is the 
mass of the host halo enclosed within the satellite's distance vm- 



As our simulation uses static dark matter (DM) halos, the structural change due 
to the interaction (e.g. mass-loss and memory effect in the perturbed region) is 
not included. As a first test we decided to directly implement mass-loss to our 
restricted JV-body code. At the time of the pericenter we are instantaneously re- 
ducing the mass of both interacting galaxies by the same factor. Certainly this is 
a very rude and heuristic approach, but it will already give a hint, if the short- 
comings of the dynamical friction equation for equal mass merger can be reduced 
by including mass-loss. 

Another important fact is, that the internal structure of the involved systems 
has a large impact on the merging process. E.g. the merging time is strongly 
correlated to the velocity dispersion. Therefore we exchanged our static isothermal 
halo by a static NEW (Navarro et al \TW7}i halo. 



Merging ability for restricted N-body 
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3 Results 
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Fig. 1. Orbital decay of self-consistent equal mass mergers with different halo properties. 
Isothermal halos with different velocity dispersions (grey lines, decreasing order: solid, 
dashed, dash-dotted) are compared with a disk merger (Kuijken & Dubinski I1995|) . set 
up as one pro- and one retrograde disk (black solid line). All models have been set on an 
initially parabolic orbit with rhaio/^'min = 5 and evolved with gyrfalcON (Dehnen l2000p 
tree-code. 

It is shown in Fig.[Tl that the internal structure of a halo has a strong impact on 
the merging time. If the velocity dispersion within an isothermal halo is reduced 
by a factor of 1.22, the merging time can drop by a factor of up to 5. If disk- 
like galaxies are used instead of pure isothermal halos, the merging time gets 
even smaller. In order to achieve reasonalbe merging times, we have decided to 
use NFW halos for future restricted N-body models, in addition they are able 
to describe the halo properties of a disk galaxy better than isothermal spheres. 
In order to quantify both processes from a stellar-dynamical point of view, we 
performed simulations with NFW profiles and the direct ad-hoc implementation 
of mass-loss described above. The results are shown in Fig.[2l for an appropriate 
strength of the frictional force, the radial decay can be modeled convincingly, 
though the implemented mass-loss model is very simple. At the moment we are 
improving the approach by measuring the mass-loss in self-consistent simulations 
and gauging the simplified approach. 

4 Conclusions and outlook 

In order to model observed major merger systems, we have implemented dynamical 
friction and a heuristic approach for the mass-loss. We have found, that NFW 
halos are more suitable to serve as host halos for our restricted JV-body models 
than isothermal halos (merging time, description of disk galaxy halo). Our results 
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Fig. 2. Remodeling of the radial decay of an equal mass disk merger (same as shown in 
Fig-Ill black solid line). Restricted models used a mass and distance dependent Coulomb 
logarithm. Isothermal halo, no mass-loss (grey dashed line); NFW halo, 1/3 of the total 
mass is lost at each pericenter (black dash-dotted line). 

show, that implementation of mass-loss is necessary when modeling 1 : 1 mergers 
in restricted Af-body models. At the moment we improve our implementation of 
mass-loss, which will be calibrated by measuring the mass-loss of self-consistent 
simulations. 

This work was supported by the German Science Foundation (DFG) under the grant TH 511/9-1, 
which is part of the DFG priority program 1177. 
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